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ABSTRACT 

Mathematical modeling and numerical solution are presented for the flow of a Micro polar fluid under the 
effect of magnetic field in an inclined channel of half width 'a' under the considerations of low Reynolds 
number. A long wavelength approximation is used to solve the flow problem. The effect of various parameters 
on the pumping characteristics is discussed through graphs. 
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I. Introduction 

Peristaltic pumping is a form of fluid transport, generally from a region of lower to higher pressure, 
by means of a progressive wave of area contraction or expansion which propagates along the length of a tube 
like structure. Peristalsis occurs naturally as a means of pumping biofluids from one place of the body to 
another. This mechanism occurs in the gastrointestinal, urinary and reproductive tracts and many other 
glandular ducts in the living body. The early reviews of Ramachandra Rao and Usha [1] , Jaffrin and 
Shapiro[2] , Man ton [3], Brasseur et al. [4], Srivastava and Srivastava [5], Provost and Schwarz [6], Shukla 
and Gupta [7], Misra and Pandey [8] ,Rao and Rajagopal [9],Kavitha et.al[10] Vajravelu et al. [11-15], Subba 
Reddy[16,17], Srinivas [18,19] deal with the peristaltic transport of viscous fluids through tubes and channels 
having impermeable flexible walls. 

Eringen [20,21] reported the theory of micropolar fluids in which the fluid micro elements undergo 
rotations without stretching. Micropolar fluids are superior to the Navier-Stokes fluids and they can sustain 
stresses and body couples. Here the micro particles in the volume Av rotate with an angular velocity about the 
centre of gravity of the volume in an average sense and is described by the micro rotation vector. The 
micropolar fluids can support stress and body couples and find their applications in a special case of fluid in 
which micro rotational motions are important. Airman and Cakmak [22] discussed three basic viscous flows of 
micropolar fluids. They are Couette and Poiseuille flows between two parallel plates and the problem of a 
rotating fluid with a free surface. The results obtained are compared with the results of the classical fluid 
mechanics. Srinivasacharya et al.[23] made a study on the peristaltic pumping of a micro polar fluid in a tube. 

Magnetohydrodynamics (MHD) is the science which deals with motion of highly conducting fluids in 
the presence of a magnetic field. The motion conducting fluid across the magnetic field generates electric 
currents which change the magnetic field, and the action of the magnetic field on these currents gives rise the 
mechanical forces which modify the flow of the fluid. (Ferraro V C A, [24]). The effect of moving magnetic 
field on blood flow was studied by Stud et al (Stud et al.[25]), and they observed that the effect of suitable 
moving magnetic field accelerates the speed of blood. 

Krishna Kumari et.al [26] studied the peristaltic pumping of a Casson fluid in an inclined channel 
under the effect of a magnetic field. Krishna Kumari et.al [27] studied the peristaltic pumping of a Jeffrey fluid 
in a porous tube. Ravi Kumar et.al [28] studied the unsteady peristaltic pumping in a finite length tube with 
permeable wall. However, the peristaltic transport of micropolar fluids in an inclined channel in the presence 
of magnetic field has not been studied. 
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In view of this, we considered the peristaltic pumping of a micropolar fluid in an inclined channel 
under the effect of magnetic field. This mathematical model may be useful to have a better understanding of 
the physiological systems such as blood vessels. 



H. MATHEMATICAL FORMULATION 

Consider the peristaltic pumping of a micropolar fluid in an inclined channel of half -width 'a'. A 
longitudinal train of progressive sinusoidal waves takes place on the upper and lower walls of the channel. For 
simplicity we restrict our discussion to the half-width of the channel as shown in Fig.l. The wall deformation 
is given by 



2k 

a + b sin ( X 

A 



ct) 



H (X ,t) 

where b is the amplitude, X is the wavelength and c is the wave speed. 
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Fig. 1 Physical Model 

Under the assumption that the channel length is an integral multiple of the wavelength A and the 
pressure difference across the ends of the channel is a constant, the flow becomes steady in the wave frame 
(x , y ) moving with velocity c away from the fixed (laboratory) frame ( X , Y ) . The transformation between 
these two frames is given by 

x = X - ct ; y = Y ;u(x, y) = U (X - ct ,Y ); v(x, y) = V (X - ct ,Y ) (2) 
where U and V are velocity components in the laboratory frame and u , v are velocity components in the 

wave frame. In many physiological situations it is proved experimentally that the Reynolds number of the flow 
is very small. So, we assume that the flow is inertia-free. Further, we assume that the wavelength is infinite. 
Using the non-dimensional quantities. 
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The non-dimensional form of equations governing the motion (dropping the bars) are 
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is coupling number, Q is the micro rotation velocity, u is the velocity, ju is the viscosity of 



the fluid, k is the micropolar viscosity, m is the micropolar parameter, p is the fluid pressure, m is the 
Hartmann number. 
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The corresponding non-dimensional boundary conditions are 



^Uo aty = (5) 

dy 

^ = at y = (6) 

dy 



Da a " at y = h(x) (7) 

a dy 

^L^Laty = h(x) (8) 



n = -l - - 

a dy 

where Da is the Darcy number and a is the slip parameter. 

III. Solution of the problem 

The general solution of (3) and (4) using the boundary conditions (5) -(8) is given by 

u = L w cosh( L 3 y)( L f L lt + L 9 )c l + —e Hy c x {L u - l)+ (e L *' + e~ L ' s )( L u + l)c, + L l5 {L w L t cosh( L 3 y 



■ ' ■[, • < 1 ,9f (9) 
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L, = L,(L,L 5 - L 3 ); L 8 = L 4 (L,(L 3 - L 4 ) - (L 3 - L 4 )); = L 4 (L 5 (L 3 + L„) + (L, - L 4 )); 

L 12 = L 8 L 10 (cosh[ L,/i] - L n L 3 sinh[ L : h]) - L 4 (( 1 / L 3 ) + L n ) exp[ L 3 /i] + (1 - L _ L 4 ) exp[ L 4 ft]; 
L 13 = L 9 L 10 (cosh[ L,A] - L U L } sinh[ L,A]) + i 4 ((1 / £ 3 ) + L u ) exp[ L 3 A ] - (1 - L L 4 ) exp[ -L 4 fc]; 

L 14 =L 13 /L 12 ; L 15 = L 3 L 10 (L 5 -L 3 )(L g L 14 - L,)(sinh[ L 3 ft ] - (L 3 L„ ) cosh[ L 3 h]); 

L i6 = L 4 L 5 L i4 ( ex P[ i 3 /i] + exp[ L 4 /i]); L„ = L 4 L, 4 (L 3 exp[ L 3 /j]+L 4 exp[ L 4 /i]); 

L i8 =- L 4 - L 5 ( ex P[ i 3 /j] + exp[ -L 4 /i]); L 19 = L 4 (L 3 exp[ L 3 /j]+L 4 exp[ -L 4 /i]); 

L 20 = L 3 L 8 L |0 (L 5 -L 3 )(L 3 L U cosh[ L 3 A] - sinh[ L,/.]); 

L 21 = L 4 L 5 (L u (L, exp[ L 3 /i] + L 4 exp[ L 4 ft ]) - (exp[ L 3 h] + exp[ L 4 ft ])); 

L, 2 = L 4 (L n (L, exp[ L 3 h] + l\ exp[ Z- 4 ft]) - (l\ exp[ L 3 ft] + l\ exp[ L 4 ft])); 

L 23 = L 4 L 5 L 11 ( L 14 ( L 3 eX P[ L 3 / '] + L 4 eX P[ L 4 /i ]) ~ ( L 3 eX P[ L 3 A ] _ L 4 eX P[ -^4^]))! 

L 24 = L 4 L n (L 14 (L 3 exp[ L 3 /j] + L 4 exp[ L 4 ft]) + (L 3 exp[ L } h] - l\ exp[ -L 4 ft])); 

L 25 = L 15 + L 16 ~ L 17 + L 18 _ L I9 _ L 23 + L 24 ' L 26 = L 20 + L 21 ~ L 22 > 

^,=(-l/i 25 i 6 )-(i 15 i 26 IL 15 ) 
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The volume flux q through each cross section in the wave frame is given by 

h 

q = | udy (11) 

o 

The pressure gradient is obtained from equation (11) 

dp I dx = (M 2 - (1 - NM q - S 6 )/(S 7 - hL 12 ))L l2 - (M 2 + rj sin 0) l{\ - N ); (12) 

(S[2)= (L lQ IL 3 )Sinh (L 3 h), S 2 = i 4 (/'' - 1)/L 2 ,, S 3 = {e L '" - 1) / l\, S 4 = 5,L g - 5 2 + S 3 , 
S 5 = Sl L 9 -S 2+ S 3 , S 6 =5 4 (1-L 14 )/L 6 L 25 , S 7 = (L 26 (S 4 (1 - L, 4 ) + S 5 ) / L 25 

The time averaged flow rate is 

Q = q + 1 (13) 

IV. THE PUMPING CHARACTERISTICS 

Integrating the equation (12) with respect to over one wave length, we get the pressure rise (drop) 
over one cycle of the wave as 

i 

, dp 

Ap = f dx (14) 

{dx 

The dimensionless frictional force F at the wall across one wavelength in the inclined channel is 

given by 

F = f h{- — ]dx (15) 
J i a i 



ii 



V. Results and discussions 

The variation of pressure rise A p with time averaged flow rate for different values of tr is shown in 

Fig. 2. It is observed that for a given Q , A p decreases as the slip parameter a increases in the pumping and 
free pumping regions. The opposite behavior is observed in co-pumping region. And also for a given Ap the 

flux Q increases with increase in a ■ 

The variation of pressure rise with time averaged flow rate for different values of Hartmann 
number M is shown in Fig. 3. It is observed that for a given Q , Ap decreases for a decreasing M in pumping 

and free pumping regions. For a given A p the flux Q depends on M and it increases with 
increasing M .The variation of pressure rise with time averaged flow rate for different values of micropolar 
parameter M is shown in Fig.4. It is observed that for a given Q, Ap increases for a increasing Tfi in pumping 

and free pumping regions. For a given A p the flux Q depends on m and it increases with increasing m . 

The effect of the inclination angle # on pumping characteristics is shown in Fig. 5. It is observed that 
for a given Q , Ap increases as the angle of inclination 6 increases. Also for a given 6 , Ap increases as Q 
increases. The variation of pressure rise with time averaged flow rate for different values of Darcy number Da 
is shown in Fig. 6. It is observed that for a given Q , Ap increases for a increasing in Da pumping and free 

pumping regions. For a given A p the flux Q depends on Da and it increases with increasing Da. 

The effect of coupling parameter on the pumping characteristics is shown in Fig. 7. We observed that 
the large the coupling number N , the pressure rises against which the pumping works. For a given Q the 

pressure difference increases with increase in N .The effect of amplitude ratio on pumping characteristics is 
shown in Fig. 8. It is observed that the large the amplitude ratio, the greater the pressure rise against which the 

pump works. For a given A p , the flux Q depends on h and it increases with increasing h . 
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The effect of rj on pumping characteristics is shown in Fig. 9. It is observed that for a given Q , Ap 

increases as rj increases. Also for a given 77 , Ap increases as Q increases. Figs. 10 to 14 are drawn to study 

the effect of various parameters on the microrotation velocity. From Fig. 10 it is observed that an increase in the 
slip parameter decreases the microrotation velocity. From Fig. 11 it is noticed that decrease in the darcy 
number decreases the microrotation velocity. From Fig. 12 it is observed that increase in M increases the 
microrotation velocity. Similarly increase in coupling parameter increases the microrotation velocity and is 
shown in figure. The effect of micropolar parameter on the microrotation velocity is shown in Fig. 13. It can be 
seen that the decrease in m decreases the microrotation velocity. 

VI. CONCLUSIONS 

Mathematical modeling of the peristaltic pumping of a Micropolar fluid under the effect of a magnetic 
field in an inclined channel is done in this paper. The following are the conclusions drawn from this. 

1. Pumping decreases as the slip parameter a increases in the pumping and free pumping regions. The 
opposite behavior is observed in co-pumping region. 

2. For a given time averaged flow rate, the pressure difference decreases for a decreasing magnetic 
parameter. 

3. The increase in micropolar parameter , increases the pumping in all the pumping regions. The same 
phenomenon is observed for the angle of inclination, Darcy number also. 

4. The effect of various parameters on the microrotation velocity is studied. An increase in the slip 
parameter decreases the microrotation velocity. A decrease in the Darcy number decreases the 
microrotation velocity. An increase in the magnetic parameter increases the microrotation velocity. 
Similarly increase in coupling parameter increases the microrotation velocity. 
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Fig. 10. Variation of £1 with y for different values of slip parameter a . 
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